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We briefly review Bethe Ansatz solutions of the integrable open spin-^ XXZ quan- 
tum spin chain derived from functional relations obeyed by the transfer matrix at 
roots of unity. 



1. Introduction 



^H. A long standing problem has been to solve the open spin-^ XXZ quantum 

(-H I spin chain with general integrable boundary terms, defined by the Hamil- 

tonian ^'^ 



Af-l 



^ : ^ = ^{ II + + (1) 

+ shr^l^ctha- th/?_crj + cscha- sech/?_(ch6'_cr^ + iskO-Oi) 

— cth a+ th P+'^N + csch a+ sech /3+ ( ch 6'+(T^r + « sh S+cr^) | , 

where tr^ , cr^ , tr^ are the standard Pauli matrices, rj is the bulk anisotropy 
parameter, a± , (3± , 9± are arbitrary boundary parameters, and N is the 
number of spins. Determining the energy eigenvalues in terms of solutions 
of a system of Bethe Ansatz equations is a fundamental problem, which 
has important applications in integrable quantum field theory as well as 
condensed matter physics and statistical mechanics and perhaps also 
string theory. (For an introduction to Bethe Ansatz, see e.g. Refs. 4, 5, 6.) 

The basic difficulty in solving (1) is that, in contrast to the special case 
of diagonal boundary terms {i.e., a± or P± ±oo, in which case Ti. has 
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a ?7(1) symmetry) which was solved long ago a simple pscudovacuum 
state does not exist. For instance, the state with all spins up is not an 
eigenstate of the Hamiltonian. Hence, many of the techniques which have 
been developed to solve integrable models cannot be applied. 

We observed some time ago that, for bulk anisotropy parameter 
values 

.= 1,2,... (2) 

(hence q = is a. root of unity, satisfying qP~^^ = — 1) and arbitrary values 

of the bomidary parameters, the model's transfer matrix t(u) (sec Sec. 2) 
obeys a functional relation of order p+1. For example, for the case p = 2, 
the functional relation is 

t{u}t{u + r])t{u + 2r]) - 6{u - i])t{u + rj) - S{u)t{u + 2r}) 

-5{u + ii)t{u) = }{u), (3) 

where 5{u) and f(u) are known scalar functions which depend on the bound- 
ary parameters. (Expressions for these functions in terms of the boundary 
parameters in (1) arc given in Ref. 18.) Similar results had been known 
earlier for closed spin chains. 

By exploiting these functional relations, we have obtained Bethe Ansatz 
solutions of the model for various special cases of the bulk and boundary 
parameters: 

(i) [Rcfs. 15, 16, 17] The bulk anisotropy parameter has values (2); 
and the boundary parameters satisfy the constraint 

a_ + /?_ + a+ + /3+ = ±{9- - 9+) + rjk , (4) 

where k S [-{N +1),N +1] is even (odd) if N is odd (even), 
respectively. 

(ii) [Ref. 18] The bulk anisotropy parameter has values (2) withp even; 
and either 

(a) at most one of the boundary parameters is nonzero, or 

(b) any two of the boundary parameters {a_, a+, /?_,/?+} are ar- 
bitrary, the remaining boundary parameters from this set are 
either t] or in/2, and 6- = 9+. 

(iii) [Ref. 19] The bulk anisotropy parameter has values (2) with p 
odd; at most two of the boundary parameters {q;_, a+, /?_, /?+} 
are nonzero, and 9- =9+. 
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All of these cases have the property that the quantity A(u), defined by 

p 

^{u) = f{uf -^^6{u + jii), (5) 

3=0 

is a perfect square. 

Solutions for generic values of the bulk anisotropy parameter and for 
boundary parameters obeying a constraint similar to (4) have been dis- 
cussed in Refs. 20, 21, 22. 

Here we briefly review our results for the cases (i) - (iii). 



2. Transfer matrix 

The transfer matrix t(u) of the open XXZ chain with general intcgrable 
boundary terms, which satisfies the fundamental commutativity property 
\t{u) , t{v)\ = 0, is given by ^ 

t{u) = tro K^{u) To{u) Kq {u) fo{u) , (6) 

where To{u) and Tb(u) are the monodromy matrices 

To{u) = Ron{u) ■ ■ ■ Roi{u) , fo{u) = Roi{u) ■ ■ ■ Ron{u) , (7) 

and tro denotes trace over the "auxiliary space" 0. The R matrix is given 
by 



R{u) = 



/sh{u + r)) \ 
shu sh?7 
sh ?7 sh M 
V sh{u + r])J 



(8) 



where r] is the bulk anisotropy parameter; and K^{u) are 2x2 matrices 
whose components are given by ^'^ 



and 



-^11 ('") = 2 (sh a- ch (3-chu + ch a_ sh /3_ sh u) 
iV'^(u) = 2 (sha_ ch/3_ chu — cha_ sh/3_ shu) 
ii'-(u) = e^- sh2u, K^iiu) = e"^- sh2u. 



K+{u)=K-{-u-r]) 



(9) 



(10) 



where ,6^ arc the boundary parameters. The Hamiltonian (1) is 

proportional to t'{0) plus a constant. 
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The transfer matrix also has in periodicity 

t{u + m) = t{u) , (11) 

crossing symmetry 

t{-u-7j)=t{u), (12) 
and the asymptotic behavior 

g«(2Ar+4)+77(Ar+2) 

i(u)~-ch(^_-^+) ^^nTi ^ + --- u^oo. (13) 

3. Case (i) 

Our main objective is to determine the eigenvalues A(m) of the transfer 
matrix t{u) (6), from which the energy eigenvalues can readily be computed. 
The functional relations for the transfer matrix (e.g., (3)) evidently imply 
corresponding relations for A(u). Following Ref. 23, we observe that the 
latter relations can be written as 

detAl(u) = 0, (14) 

where M{u) is the {p+1) x (p + 1) matrix 

/ A(n) ... -h{u) \ 



Miu) 



(15) 



V ...-h{u + inj)A{u+pr,)J 



h{u) 

if there exists an ZTr-periodic function h{u) such that 

/(.) = nM» + ;,) + rig±g^ (le) 

To solve for h{u), we set z{u) = 11^=0 '^(^ + jv)^ and observe that (16) 
implies that z{u) satisfies a quadratic equation 

p 

z{uf - z{u)f{u) + W 5{u + j7?) = , (17) 

whose solution is evidently given by 

z{u) = \ {j{u) ± , (18) 
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where A(m) is defined in (5). If the boundary parameters satisfy the con- 
straint (4), then A(u) is a perfect square, and two solutions of (16) are 

sh(2u + 2r]) 



4sh^'^ {u + r])- 



sh(2w + rj) 

X sh{u ± a-) ch(u ± /3_) sh{u ± a+) ch(u ± /3+) . (19) 

Let us now label the corresponding matrices (15) by A^^^'(u). 

The condition (14) implies that Al^^^(u) has a null eigenvector v^'^^u), 



M^^\u) t;W(w) = 0, 
Note that the matrix M.^'^^u) has the symmetry 

5A^W(u)S'-^ = Al(=^)(w + r?) . 



(20) 
(21) 



where 



/O 1 . . 
1.. 



S = 



0\ 




... 1 
\l ... 0/ 



SP+'^ = I . 



(22) 



It follows that the null eigenvector v^'^\u) satisfies Sv^'^^u) = v'-'^^u + r]). 
Thus, its components can be expressed in terms of a function 

(u) = (0(±) (u) , Q(±) (« + 77) , . . . , (u + pr,)) , (23) 

with Q^^'>{u + iir) = Q'^^^u). We make the Ansatz 

M<±) 

Q(±) (u) = JJ sh(w - «f ^) sh(w + «f ^ + r?) , (24) 

which has the crossing symmetry Q^^\u) = Q^^''(— u — 77). Substituting 
the expressions for A4'^^\u) (15) and v^^\u) (23) into the null eigenvector 
equation (20) yields the result for the transfer matrix eigenvalues 

The asymptotic behavior (13) implies that M'-'^^ = ^{N — l±k), where k is 
the integer appearing in the constraint (4). Analyticity of the eigenvalues 
(25) implies the Bethe Ansatz equations 



/i(±)(-Mf'-r7) " Q(±)(wf' -r?) ' 



i = l,...,MW. (26) 
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In short, for case (i), the eigenvalues of the transfer matrix (6) are given by 
(25), where h^^\u) and Q^^'^Hu) are given by (19), (24) and (26). 

In Ref. 17, we have verified numerically that this solution holds also 
for generic values of r], which is consistent with Refs. 20, 21, 22; and that 
this solution gives the complete set of 2^ eigenvalues. To illustrate how 
completeness is achieved, let us consider the case = 4. The integer k 
in the constraint (4) must therefore be odd, with — 5 < < 5. The six 
possibilities are summarized in Table 1. 



Tabic 1. Completeness for Af = 4. For each k, there are 2^ eigenvalues. 



k 


# eigenvalues given by A(+'(u) 


# eigenvalues given by A' ^(u) 


5 


16 





3 


15 


1 


1 


11 


5 


-1 


5 


11 


-3 


1 


15 


-5 





16 



4. Case(ii) 

A key feature of case (i) is that the quantity A{u) (5) is a perfect square. 
We therefore look for additional such cases. For p even, we find that A(u) is 
also a perfect square if either (a) at most one of the boundary parameters 
is nonzero; or (b) any two of the boundary parameters {a-,a+, P+} 
are arbitrary, the remaining boundary parameters from this set are either 
r] or i7r/2, and 0^ = 0-^-. For definiteness, we foctis here on the stibcase (b) 
with a± arbitrary, P± = r] and A'' even. Unfortunately, the resulting z{u) 
(18) is not consistent. To surmount this difficulty, we use a matrix M{u) 
which is different from (15), namely 

M{u) = (27) 

/ A(w) -h{u) ... -h{-u+pr])\ 

-h{-u) K{u + prj) -h{u + prf) . . . 

X-h^a + p^rj) ... -h{-u - p{p - 1)7]) A{u + p'^r]) / 

where h{u) is 2j7r-periodic. This matrix has the symmetry 



SM{u)S-^ = M{u + pr]), 



(28) 
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where S is given by (22). By arguments similar to those used in Sec. 3, we 
find that the transfer matrix eigenvalues are given by 

A{u) = h{u) +h{-u + pr]) , 29 

Q{uj Q[u) 

where h(u) is given by 

h{u) = 4sh-(. + oh\u - n) (30) 



X sh(M — a_) sh('u + a+) 



ch (i(u + a- + 77)) ch (i(u - a+ + 77)) 



ch (^(u - a_ - T])) ch (^(w + a+ — r])) ' 
and Q{u) is given by 

= Ilsh Q(m-Mj)^ sh Q(u + Wj -pr?)^ , (31) 

with M = A'' + 2p + 1; and the Bethe Ansatz equations are 



h{uj) _ Q{uj-pri) 



j^l,...,M. (32) 



h{-Uj +pri) Q{uj +pr]) ' 

We have verified numerically the completeness of this solution. The other 
subcases (a) and (b) are mostly similar. ^ 



5. Case(iii) 

For p odd, we find that the quantity A(u) (5) is also a perfect square if at 

most two of the boundary parameters {a_, a+, /?+} are nonzero, and 
6- — For definiteness, we focus here on the case with a± arbitrary, 
/3-I- = and even. As in case (ii), the resulting z{u) (18) is not consistent. 
To surmount this difiiculty, we again use a matrix M{u) which is different 
from (15), namely 

M{u) = 
( Hu) 

\-/i(2)(u - ?y) ... -/i(i)(u+ (p- 1)77) A(M + pr/) y 



^The exception is tiic subcase (a) witli 8± nonzero, for which Q{u) = llj=i sh(u - Uj), 
which is not crossing symmetric. See Sec. 3.3 in Ref. 18. 



6{u—rj) 
■/((2)(u-r)) 





(33) 
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where h^^^ (u) and h^"^^ (u) are i7r-periodic. It has the reduced symmetry 

TM (w)T-i =M{u + 2r]), (34) 

where T = S"^, and S is given by (22). (While (21) implies (34), the 
converse is not true.) The condition det Al('u) = implies that Ad{u) has 
a null eigenvector v{u), 

Miu)v{u)=0, (35) 

where v{u) satisfies Tv{u) = v{u + 2t]). Thus, its components are expressed 
in terms of two independent functions Qi{u), Q2{u): 

v{u) = {Qi{u),Q2{u),... ,Qiiu-2r]),Q2{u-2rj)) . (36) 

We make the Ansatze 

Ml 

Qi{u) = Yl sh(u - 4'') sh(w + + ri) , 

3=1 

M2 

Q2 (u) = [| sh(w - uf ) sh{u + uf + 2,ri) . (37) 

Substituting the expressions for M.{u) (33) and v{u) (36) into the null 
eigenvector equation (35) yields two expressions for the transfer matrix 
eigenvalues, 

5{u) Q2{u) S{u-v) Q2{u-2n) 
^ ' hW{u)Qi{u) ^ h(^){u-r]) Qiin) ' 

= h^'\u-,)%^^+h^^\u)9l^ (38) 
Q2{u-r)) Q2{u-r]) 

Analyticity of these expressions leads to the Bethe Ansatz equations 
5{uf^) /i(2)(«W -rj) _ Q2{uf^ - 2r]) 



(1) , i = 1,2,... ,Mi, 

i = l,2,... ,M2. (39) 



5{uf -^)h^^){uf) Q2{U 



We expect that there are sufficiently many equations to determine all the 

zeros {u^j^\ u'"j^^ } of Qi{u) ,Q2{u), respectively. Functions h^^\u) (with 
h'^^^u) = h'^^\—u — 2r])) which ensure the condition det M.{u) = are 
given by 

/i(i)(u) = 4sh2^(w + 2?7), M2 = iiV + i(3p- 1), Mi = Ma + 2 , 
p = 3,7,ll,... (40) 
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and 



-2 ch(2u) sh^ u sh^^{u + 2r)) , Mi 
p = 9,17,25,... 

Ml 

p = 5,13,21,... 
2 ch(2u) sh^ u sh^^ (u + 2ry) , Mi 



2ch(2u) sh^ush^^{u + 2r]) , 



M2 = ^N + 2p-l, 

M2 = 

M2 = \N + 2, 



p=l. 



We have verified numerically the completeness of this solution. Similar 
results hold for the case with a_, /?_ arbitrary and a+ = /?+ = 0, etc. 

We observe that this solution represents a generalization of the famous 
Baxter T — Q relation ^, which schematically has the form 



t{u) Q{u) = Q{u') + Q{u") . 

Indeed, our result (38) has the structure 

tiu) Qiiu) = Q2{u') + Q2{u") , 
t{u) Q2(u) = Qi{u') + Qi{u") . 



(42) 



(43) 



Such generalized T — Q relations, involving two or more independent Q('u)'s, 
may also appear in other integrable models. 



6. Conclusions 

We have seen that Bethe Ansatz solutions of the open spin-| XXZ quantum 

spin chain are available for the cases (i)-(iii), for which the quantity A(?i) 
(5) is a perfect square. There may be further special cases for which A(m) 
is a perfect square, in which case it should not be difficult to find the 
corresponding Bethe Ansatz solution. Our solution for case (iii) involves 
more than one Q{u). This is a novel structure, which should be further 
understood. The general case that A(u) is not a perfect square and/or 
that ?7 7^ i-K/{p + 1) also remains to be understood. 
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